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We  studied  A1B  nanowires  as  hydrogen  storage  materials  based  on  density  functional  theory  and  Rigged 
QED  theory.  In  this  paper,  we  focused  on  the  adsorption  energy  and  the  electronic  structure  of  models.  A1B 
nanowire  models  are  compared  with  an  A1  nanowire  model  and  A1B2  crystal  structure  in  terms  of  density 
of  states,  electron  density,  kinetic  energy  density,  tension  density  and  stress  tensor  density.  These  results 
revealed  A1B  nanowires  do  not  have  the  conductivity,  while  the  A1  nanowire  and  A1B2  bulk  have  it.  It  was 
also  shown  that  the  stabilization  energies  of  A1B  nanowires  for  the  hydrogen  adsorption  are  larger  than 
that  of  A1  nanowire.  Adsorped  hydrogens  are  more  stable  in  the  A1B  nanowires  than  the  A1  nanowire. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nowadays,  hydrogen  energy  is  focused  on  as  a  clean  energy 
resource,  taking  the  place  of  fossil  fuel.  However,  hydrogen  is  gas 
in  room  temperature  and  1  atm,  and  thus  hydrogen  is  difficult  to 
store  safely.  One  of  methods  of  hydrogen  storage  is  using  hydro¬ 
gen  adsorption.  For  the  hydrogen  adsorption,  higher  weight  percent 
storage  is  required  for  hydrogen  to  take  rank  with  fossil  fuel.  How¬ 
ever,  high-weight-percent  hydrogen  storage  has  not  been  achieved 
for  most  of  hydrogen  storage  materials.  To  overcome  the  difficulty, 
nanostructures  are  focused  on.  Nanostructures  have  larger  surfaces 
than  usual  bulk  structures.  In  this  point,  they  have  a  possibility 
to  store  much  hydrogen.  Thus  many  studies  are  reported  about 
hydrogen  storage  nanostructures  [1-3]. 

In  nanomaterials,  we  particularly  focused  on  one-dimensional 
structures  such  as  a  carbon  nanotube  (CNT)  and  A1  nanowire 
[1,4,5].  Makita  et  al.  [4]  calculated  optimized  geometries  of  alu¬ 
minum  nanowires  based  on  gold  nanowires  using  first-principles 
electronic  structure  calculations.  Kawakami  et  al.  [5]  showed  that 
hydrogen  is  adsorbed  on  a  pentagonal  Al  nanowire  model  by  the 
separation  of  a  H2  molecule  into  two  H  atoms.  Nakano  et  al.  [1] 
suggested  to  wrap  Al  species  in  carbon  materials  to  enhance  the 
hydrogen  adsorption  on  CNT.  Note  that  Suzuki  et  al.  [6]  reported 
that  Al  whiskers  are  actually  fabricated  by  glancing  angle  deposi¬ 
tion  on  a  high  temperature  (HT-GLAD)  substrate. 

For  a  hydrogen  molecule  adsorption  on  these  hydrogen  storage 
materials,  the  low  adsorption  and  dissociation  energy  under  1.0  eV 
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is  for  the  usage  as  a  fuel  cell.  However,  these  energy  are  generically 
larger  than  1.0  eV  for  these  nanostructures.  As  a  solution  strategy  to 
overcome  the  high  adsorption  and  desorption  energy,  we  suggest 
the  use  of  the  electronic  and  magnetic  force.  Thus  in  this  paper, 
we  focus  on  hydrogen  adsorption  abilities  and  the  electronic  struc¬ 
tures  of  the  A1B  and  Al  nanowires.  A1B  nanowires  are  expected  to 
be  a  better  candidate  for  hydrogen  storage  because  of  light  weight 
of  B  comparing  to  Al.  We  carry  out  first-principles  calculations  to 
optimized  geometries,  analyze  electronic  structures  and  clarify  the 
influence  of  substitution  of  Al  with  B  for  the  electronic  structures. 
In  addition,  we  analyze  details  of  the  electronic  structures  in  terms 
of  the  quantum  mechanical  energy  densities  based  on  the  regional 
density  functional  theory  [7-13]. 

2.  Computational  methods 

2.1.  Calculation  model  and  method 

Two  A1B  nanowire  models  for  our  calculation  are  shown  in 
Figs.  1  and  2.  Fig.  1(a)  shows  A1B2  bulk  models  [14].  Fig.  1(b)  and  (c) 
shows  the  hexagonal  nanowire  models  based  on  A1B2  bulk  mod¬ 
els.  Fig.  1(d)  and  (e)  are  A1B2  bulk  surface  models,  which  are  called 
as  the  surface  model  1  and  model  2,  respectively.  These  models  are 
also  based  on  A1B2  bulk  models.  In  this  paper,  the  small  nanowire  of 
Fig.  1  (c)  is  considered  to  be  a  hydrogen  storage  material.  Other  mod¬ 
els  are  used  to  compare  with  nanowire  models.  In  Fig.  2,  another  A1B 
nanowire  and  Al  nanowire  models  are  shown.  This  A1B  nanowire 
model  is  composed  of  pentagonal  B  rings  and  an  Al  monoatomic 
wire.  Fig.  2(a)  shows  the  Al  nanowire  whose  geometry  was  cal¬ 
culated  by  Makita  et  al.  [4]  consisting  of  pentagonal  rings  and 
monoatomic  wire.  Pentagonal  B  rings  are  placed  with  rotating  by 
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Fig.  1.  (a)  A1E>2  bulk  model,  (b)  A1B  larger  hexagonal  nanowire  model,  (c)  A1B  hexagonal  nanowire  model,  (d)  A1B2  bulk  surface  model  1,  and  (e)  A1B2  bulk  surface  model  2. 


36°  for  each  ring.  The  A1B  pentagonal  nanowire  shown  in  Fig.  2(b) 
is  based  on  this  Al  nanowire.  This  A1B  nanowire  is  also  considered 
to  be  a  hydrogen  storage  material. 

The  boundary  condition  of  these  models  is  imposed  as  peri¬ 
odic  one.  In  our  calculation,  two  kinds  of  supercell  are  used. 
8.36  A  x  8.36  A  x  4c  A  supercell  is  considered  for  A1B  nanowires  and 
10.6  A  x  10.6  A  x  8c  A  supercell  is  considered  for  Al  nanowire.  The 
parameter  c  is  the  distance  from  atoms  on  the  axis  to  the  pen¬ 
tagonal  or  hexagonal  ring.  If  the  smaller  unit  cell  is  used  for  Al 
nanowire,  the  interaction  with  next  unit  cell  is  not  negligible.  The 
supercells  for  A1B2  surface  models  shown  in  Fig.  1(d)  and  (e)  are 
5.19  A  x  6.54  A  x  10.4  A  and  5.99  A  x  6.54  A  x  10.4  A,  respectively. 

To  investigate  the  electronic  structures  of  the  systems,  vari¬ 
ational  calculations  are  carried  out  based  on  density  functional 
theory,  which  is  implemented  by  FFII98md  program  package 
[15].  Wave  functions  of  valence  electrons  are  expanded  by 
Bloch  functions,  while  atomic  cores  and  inner-shell  electrons  are 
replaced  by  Hamman’s  norm  conserving  pseudopotentials  [16]. 
The  generalized-gradient-approximation  method  is  employed  for 
exchange-correlation  interaction  formulated  by  Perdew  and  Wang 
(PW91)  [17].  The  values  of  cutoff  energy  for  wave  functions  with 


plane-wave  expansion  are  chosen  as  1088  eV  for  the  A1B  nanowires 
and  408  eV  for  the  Al  nanowire.  The  higher  cutoff  energy,  1088  eV,  is 
required  to  derive  an  accurate  result  for  B  atoms.  The  cutoff  energy 
for  bulk  surface  model  is  the  same  as  the  lower  cutoff  energy.  This 
is  because  the  vacuum  region  is  only  one  direction.  Thus  the  lower 
cutoff  energy  is  sufficient. 

The  density  of  states  (DOS)  are  calculated  using  a  density  func¬ 
tional  method  with  ADF  program  package  [18].  Electron  wave 
functions  are  expanded  with  linear  combinations  of  atomic  orbital 
basis  set.  In  our  calculation,  hybridization  of  Slater  type  double-zeta 
basis  and  numerical  one  are  employed. 

The  minimum  point  in  a  potential  energy  surface  corresponds 
optimized  structure.  In  the  optimization  calculation,  two  geometric 
parameters,  the  radius  of  nanowire  r\  and  the  length  of  the  unit  c, 
are  determined.  The  geometries  of  our  models  are  optimized  with 
keeping  their  symmetries. 

2.2.  Quantum  energy  density 

Quantum  energy  density  [7-10]  provides  a  number  of  informa¬ 
tion  about  electronic  properties.  First,  we  can  derive  kinetic  energy 
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Fig.  2.  (a)  Al  pentagonal  nanowire  model  and  (b)  A1B  larger  hexagonal  nanowire 
model. 


density  nr(r)  from  electron  wave  functions  as  follows: 


mo = 

i 


(1) 


where  vt  is  the  occupation  number  of  the  ith  state  T/q.  This  energy 
density  represents  whether  electrons  can  behave  as  classical  par¬ 
ticle.  Electrons  can  freely  move  as  obeying  classical  mechanics  in 
positive  kinetic  energy  density  region,  which  is  called  as  electronic 
drop  region  ftD.  On  the  other  hand,  in  the  negative  kinetic  energy 
density  region,  which  is  called  as  electronic  atmosphere  region  RA, 
electrons  cannot  move  freely  and  exist  by  only  quantum  tunnel¬ 
ing  effects.  The  boundary  of  these  two  regions  is  called  electronic 
interface  S. 

Stress  tensor  density  r  (r)  is  denoted  as  a  3  x  3  matrix  whose 
element  r skl(f)  is  given  by 


rskl(r)  = 


Mh(r) 

dxkdxl 


ty*(r)  djsj(r) 
dxk  dxl 


dxkdxl  1  dx i  dxk 


(2) 


The  positive  and  negative  eigenvalues  of  r  (r )  represent  the  tensile 
and  compressive  stress,  respectively.  This  tensor  clarifies  proper¬ 
ties  of  chemical  bonds,  i.e.,  covalent  bonds  and  metallic  ones  are 
distinguished.  In  this  paper,  we  show  the  largest  eigenvalues  and 
their  eigenvectors  as  the  principal  stress. 


Fig.  3.  Total  energy  curves,  (a)  A1B  pentagonal  nanowire  model  as  a  function  of  ri,  (b)  A1B  pentagonal  nanowire  model  as  a  function  of  c,  (c)  A1B  hexagonal  nanowire  model 
as  a  function  of  n ,  and  (d)  A1B  hexagonal  nanowire  model  as  a  function  of  c. 


Fig.  4.  Total  A£/Ac  as  a  function  of  Ac.  (a)  A1B  pentagonal  nanowire  model  and  (b)  A1B  hexagonal  nanowire  model. 
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Fig.  7.  ( Continued ). 


Tension  density  rs(r)  =  (rsl(r),  rS2(r),  rS3(r))  is  given  as  the 
divergence  of  the  stress  tensor  density,  r  (r): 


rsk(r)  =  - — 


h2  v- 


4mZ-^  1 

i 

dAifs*(r) 


'  Vr,C?)  9^- (?)...,. 

w— - irA*(r) 


+ 


dxk 


dxk 


(3) 


3.  Results  and  discussion 

3.1.  Geometry  ofAlB  nanowires 

Fig.  3  shows  the  potential  energy  curves  of  both  A1B  nanowires 
for  hydrogen  storage.  The  energy  curves  are  shown  depending  on 
the  parameter  r\  and  c.  The  optimized  parameter  corresponds  to 
the  minimum  energy.  The  calculation  is  done  in  detail  around 
the  stable  geometry.  In  the  case  of  the  pentagonal  A1B  nanowire, 
ri  and  c  are  determined  at  1.31  A  and  1.80  A,  respectively.  In 
this  structure,  the  bond  length  between  B  atoms  on  the  same 
ring  is  1.54  A  and  that  between  B  and  Al  atoms  is  2.23  A.  In  the 
hexagonal  A1B  nanowire,  r\  and  c  are  1.54  A  and  1.68  A,  respec¬ 


tively,  and  the  length  between  B  (B  and  Al)  atoms  is  1.54  A 
(2.28  A).  The  optimized  lengths  of  B-B  bond  are  the  same  for 
both  models,  while  other  parameters  are  different  from  each 
other. 

On  the  other  hand,  in  larger  hexagonal  nanowire  models,  the 
B-B  bond  length  is  1.67  A  and  the  B-Al  length  is  2.36  A.  Moreover 
the  interatomic  distances  in  A1B2  crystal  are  known  as  1.73  A  for 
B-B  and  2.38  A  for  B-Al,  these  bond  lengths  are  shorter  in  nanowire 
structures.  In  Al  nanowire,  the  wire  radius  r\  and  the  length  of  the 
unit  c  are  2.48  A  and  1.34  A,  respectively. 

Fig.  4  shows  AE/Ac  as  a  function  of  Ac,  where  Ac  is  the  dif¬ 
ference  between  c  and  that  of  the  stable  geometry,  and  A E  is  the 
difference  in  the  potential  energy.  In  both  simple  A1B  nanowire 
models,  the  line  can  be  approximately  considered  as  straight.  The 
gradients  of  the  lines  mean  the  stability  of  nanowires.  This  figure 
shows  that  larger  energy  is  needed  for  larger  Ac.  Therefore,  the  gra¬ 
dients  are  regarded  as  Young’s  modulus.  Thus  these  figures  indicate 
the  stabilization  of  simple  A1B  nanowires  in  the  axial  direction.  But 
the  geometry  which  is  calculated  with  keeping  symmetry  may  be  a 
meta-stable  one.  The  stability  of  nanowire  model  is  investigated  by 
the  deformation  along  the  normal  vibrational  modes,  where  the 
norm  of  the  eigenvectors  is  set  to  0.1  A.  The  normal  vibrational 
modes  are  calculated  with  Gaussian  03  in  the  molecular  system 
[19].  Atoms  of  models  are  moved  to  the  direction  of  the  imaginary 


66 


A.  Fukushima  et  al.  /  Journal  of  Power  Sources  184  (2008)  60-76 


(-7.00,-7.00,0.00)  [bohr] 


(7.00-7.00,0.00) 


(-7.00,0.00,-7.00) 


[bohr] 


(7.00,0.00,-7.00) 


(-7.00,0.00,-7.00)  [bohr] 


(7.00,0.00,-7.00) 


(-7.00,-7.00,0.00) 


[bohr] 


(e) 


(f) 


(-7.00,0.00,7.00) 


(-7.00,0.00,7.00) 


(7.00,-7.00,0.00) 


(-7.00,0.00,-7.00)  [bohr]  (7.00,0.00,-7.00)  (-7.00,0.00,-7.00)  [bohr]  (7.00,0.00,-7.00) 


Fig.  8.  Kinetic  energy  density  and  tension  density  of  nanowire  models.  The  region  in  the  circle  is  pseudopotential.  The  tension  density  is  showed  only  on  the  region  where 
the  kinetic  energy  is  positive.  Kinetic  energy  density  and  tension  density  of  the  sum  of  all  occupied  states  in  (a  and  b)  A1B  pentagonal,  (d  and  e)  A1B  hexagonal,  and  (g  and  h) 
Al  pentagonal  nanowire  models.  Light  (dark)  gray  region  means  positive  (negative)  kinetic  energy  density.  Kinetic  energy  density  and  tension  density  of  the  highest  occupied 
states  (c)  A1B  pentagonal,  (f)  A1B  hexagonal,  and  (i)  Al  pentagonal  nanowire  models. 
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Fig.  8.  ( Continued ). 


normal  vibrational  modes  and  the  total  energy  are  calculated  with 
this  models  in  the  periodic  system.  As  a  result,  these  nanowires  are 
stable  enough. 

3.1.1.  Electronic  structure  of  simple  nanowires 

First,  we  discuss  the  electronic  structure  of  nanowire  models. 
The  DOS  and  band  structure  of  simple  nanowire  models  are  shown 
in  Fig.  5.  The  Fermi  levels  are  shifted  to  the  origin  in  the  verti¬ 
cal  axes.  In  addition,  the  DOS  of  the  A1B2  bulk,  surface  and  larger 
nanowire  model  are  shown  in  Fig.  6.  Moreover  the  band  structure  of 
larger  nanowire  models  is  shown  in  Fig.  6.  The  Fermi  levels  are  also 
shifted  as  figures  of  simple  nanowire  models.  In  the  pentagonal  A1B 
nanowire,  the  band  gap  on  the  Fermi  level  is  seen  in  Fig.  5(a).  Flence 
this  indicates  that  A1B  pentagonal  nanowire  has  different  conduc¬ 
tive  property  from  A1B2  bulk  structure.  On  the  other  hands,  DOS 
and  band  structure  of  A1B  hexagonal  nanowire  in  Fig.  5(b)  is  metal¬ 
lic.  Moreover  the  DOS  and  band  structure  of  A1B2  surface  and  larger 
hexagonal  nanowire  are  shown  in  Fig.  6(a)-(c).  In  Al  pentagonal 
nanowire  model,  DOS  and  band  structure  is  metallic. 

Next,  to  study  electronic  state  of  the  simple  nanowire  models 
in  detail,  the  electron  density,  the  kinetic  energy  density,  tension 
density  and  stress  tensor  density  are  calculated.  The  electron  den¬ 
sity  on  the  xz-plane  and  xy-plane  are  shown  in  Fig.  7.  Fig.  7(a),  (b), 


(d),  (e),  (g)  and  (h)  are  the  sum  of  the  occupied  states  and  Fig.  7(c), 
(f)  and  (i)  are  the  highest  occupied  state.  Fig.  7(a),  (d)  and  (g)  show 
electron  density  on  xy-plane  and  the  other  figures  show  that  of  xz- 
plane.  In  Fig.  8,  the  kinetic  energy  density  and  the  tension  density  is 
shown.  Fig.  8(a),  (b),  (d),  (e),  (g)  and  (h)  are  the  sum  of  the  occupied 
states  and  Fig.  8(c),  (f)  and  (i)  are  the  highest  occupied  state.  The 
tension  density  is  normalized.  Fig.  8(a),  (d)  and  (g)  show  the  kinetic 
energy  density  and  tension  density  on  xy-plane  and  the  other  fig¬ 
ures  show  those  of  xz-plane.  Fig.  9  shows  the  largest  eigenvalue  and 
eigenvector  of  stress  tensor  density.  The  eigenvector  is  normalized. 
Fig.  9(a),  (c)  and  (e)  show  stress  tensor  density  on  xy-plane  and  the 
other  figures  show  that  of  xz-plane. 

First,  we  discuss  the  pentagonal  A1B  nanowire.  In  Fig.  7(a)  and 
(b),  the  electron  density  between  B  atoms  is  high.  It  shows  that 
B  atoms  are  bonded  each  other.  It  is  seen  from  Fig.  7(c)  that  the 
electron  density  of  the  highest  occupied  state  concentrates  on  B 
atoms  rather  than  on  Al  atoms.  The  highest  occupied  state  B(2s2p) 
appears  on  the  direction  perpendicular  to  the  wire  axis  like  dan¬ 
gling  bonds.  Flence,  it  is  seen  that  it  makes  bonds  not  along  the 
wire  axis  but  between  B  atoms  on  the  five-membered  ring.  There¬ 
fore,  in  the  pentagonal  A1B  nanowire,  electrons  cannot  move  along 
the  nanowire  axis  like  free  electrons.  Moreover,  it  can  be  speculated 
that  this  concentrated  electron  density  affects  H  atom  adsorption. 


68 


A.  Fukushima  et  al.  /  Journal  of  Power  Sources  184  (2008)  60-76 


(a) 


(b) 


(c) 

(-7.00,7.00,0.00) 


(d) 

(-7.00,0.00,7.00) 


O 

.Q 


l  0  005 

0.0 

I  -0.005 


(-7.00,0.00,-7.00)  [bohr] 


(7.00,0.00,-7.00) 


(-10.0,0.00,-10.0)  [bohf] 


(10.0,0.00,-10.0) 


Fig.  9.  The  largest  eigenvalues  of  stress  tensor  density  of  nanowire  models.  Red  (blue)  region  means  positive  (negative)  stress  density.  The  eigenvectors  are  normalized. 
The  largest  eigenvalues  and  eigenvectors  of  stress  tensor  density  in  (a  and  b)  A1B  pentagonal,  (c  and  d)  A1B  hexagonal,  and  (e  and  f)  Al  pentagonal  nanowire  models.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Fig.  10.  H  atom  adsorption  site.  Site  1 :  the  top  of  B(A1)  ring.  Site  2:  the  same  distance 
from  B(A1)  atoms  of  the  same  ring.  Site  3:  the  same  distance  from  B(A1)  atoms  of  the 
different  rings. 


This  tendency  is  also  seen  in  the  kinetic  energy  density  and  tension 
density.  In  Fig.  8(a)  and  (b),  the  kinetic  energy  density  concentrate 
between  B  atoms  and  in  Fig.  8(c)  the  kinetic  energy  density  of  the 
highest  occupied  state  concentrate  on  the  direction  perpendicu¬ 
lar  to  the  wire  axis  like  dangling  bonds.  In  Fig.  9(a)  and  (b),  the 
largest  eigenvalue  and  the  eigenvector  are  shown.  The  red  (blue) 


Fig.  11.  H  absorbed  simple  nanowire  models  and  A1B2  surfaces,  (a)  A1B  pentagonal 
nanowire  model,  (b)  A1B  hexagonal  nanowire  model,  (c)  Al  pentagonal  nanowire 
model,  (d)  A1B2  surface  model  1,  and  (e)  A1B2  surface  model  2. 


Fig.  12.  DOS  and  band  structure  of  nanowire  models  after  H  adsorption,  (a)  A1B  pen¬ 
tagonal  nanowire  model,  (b)  A1B  hexagonal  nanowire  model,  and  (c)  Al  pentagonal 
nanowire  model. 
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Fig.  13.  DOS  of  surface  models  after  H  adsorption,  (a)  A1B2  surface  model  1  and  (b)  A1B2  surface  model  2. 


region  means  the  positive  (negative)  eigenvalue  of  the  stress  tensor 
density.  The  positive  stress  is  seen  between  B  atoms  and  spindle 
structure  is  seen  in  the  same  region.  This  means  B-B  bond  has  the 
covalent  property. 

Properties  of  the  hexagonal  A1B  nanowire  are  similar  to  those 
of  the  pentagonal  one  as  seen  in  Figs.  7-9.  However,  in  the 
hexagonal  one,  the  highest  occupied  state  B(2s2p)  appears  not 
only  on  the  direction  perpendicular  to  the  wire  axis.  It  appears 
around  B  atoms,  which  is  in  contrast  with  that  of  the  pen¬ 


tagonal  one.  The  interaction  of  these  atomic  orbitals  is  weak, 
and  hence  narrow  band  appears  near  the  Fermi  level  as  seen 
in  Fig.  5(b).  However,  electrons  on  the  narrow  band  are  not 
conductive,  since  electrons  concentrates  on  B(2s2p)  hybridized 
orbitals  and  heavy  effective  masses  are  understood  from  the 
DOS  and  band  structure.  Thus,  the  motion  of  electrons  is 
also  constrained  only  on  the  atoms,  and  hence  this  nanowire 
is  expected  not  to  have  the  conductivity  as  the  pentagonal 
one. 
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Fig.  14.  Electron  density  of  nanowire  models  after  H  adsorption.  The  region  in  the  circle  is  pseudopotential.  Electron  density  of  the  sum  of  all  occupied  states  in  (a  and  b)  A1B 
pentagonal,  (d  and  e)  A1B  hexagonal,  and  (g  and  h)  Al  pentagonal  nanowire  models.  Electron  density  of  the  highest  occupied  states  (c)  A1B  pentagonal,  (f)  A1B  hexagonal,  and 
(i)  Al  pentagonal  nanowire  models. 
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(-7.00,-7.00,0.00)  [bohr]  (7.00,-7.00,0.00)  (-7.00,0.00,-7.00)  [bohr]  (7.00,0.00,-7.00) 


(-7.00,0.00,-7.00)  [bohr]  (7.00,0.00,-7.00)  (-7.00,-7.00,0.00)  [bohr]  (7.00,-7.00,0.00) 


(-7.00,0.00,-7.00)  [bohr]  (7.00,0.00,-7.00)  (-7.00,0.00,-7.00)  [bohr]  (7.00,0.00,-7.00) 


Fig.  15.  Kinetic  energy  density  and  tension  density  of  nanowire  models  after  H  adsorption.  The  region  in  the  circle  is  pseudopotential.  The  tension  density  is  showed  only 
on  the  region  where  the  kinetic  energy  is  positive.  Kinetic  energy  density  and  tension  density  of  the  sum  of  all  occupied  states  in  (a)  and  (b)  A1B  pentagonal,  (d)  and  (e)  A1B 
hexagonal,  and  (g)  and  (h)  Al  pentagonal  nanowire  models.  Light  (dark)  gray  region  means  positive  (negative)  kinetic  energy  density.  Kinetic  energy  density  and  tension 
density  of  the  highest  occupied  states  (c)  A1B  pentagonal,  (f)  A1B  hexagonal,  and  (i)  Al  pentagonal  nanowire  models. 


A.  Fukushima  et  al.  /  Journal  of  Power  Sources  184  (2008)  60-76 


73 


(9) 


(h) 


(-7.07,7.07,0.00) 


(-7.07,-7.07,0.00)  [bohf]  (7.07,7.07,0.00) 


(-7.07,-7.07,10.0) 


(-7.07,-7.07,-10.0) 


[bohr]  (7.07,7.07,-10.0) 


(i) 

(-7.07,-7.07,10.0) 


•  '*Th  *  _ V  r\r  "  _ 

[bohr]  (7.07,7.07,-10.0) 


Fig.  15.  ( Continued ). 


Next,  we  consider  the  pentagonal  Al  nanowire.  The  band 
structure  shows  metallic  properties  as  shown  in  Fig.  5(c).  In 
the  Al  nanowire,  the  electron  density  and  the  positive  kinetic 
energy  density  region  is  observed  in  the  whole  Al  nanowire 
(Figs.  7  and  8(g)— (i)).  Hence,  it  is  expected  that  electrons  can  be 
accelerated  by  the  external  fields  and  escape  from  the  constraint 
of  atomic  nuclei.  Therefore,  it  can  be  concluded  that  Al  nanowire 
has  conductivity.  The  largest  eigenvalue  is  almost  negative  on  the 
xz-plane  and  xy-plane.  This  shows  metallic  bond  between  Al  atoms. 
In  Al  nanowire  models,  electrons  of  HOMO  state  are  not  localized. 
Thus  it  can  be  speculated  that  A1B  nanowire  has  different  hydrogen 
adsorption  property. 


3.2.  Hydrogen  adsorption  on  AIB  nanowires 

3.2.1.  Hydrogen  adsorption  site  and  energy 

In  this  section,  hydrogen  adsorption  on  the  nanowires  is  dis¬ 
cussed.  Total  energy  is  calculated  to  find  out  the  stable  adsorption 
point  of  a  H  atom  on  the  nanowires.  Three  adsorption  points,  sites 
1,  2,  and  3,  as  shown  in  Fig.  10  are  considered.  First,  the  geometries 
of  the  nanowires  are  fixed  for  the  calculation  of  total  energy.  The 


Table  1 

The  stabilization  energy  of  H  atom  adsorption 


stabilization  energy  A E  is  defined  as 

A E  =  £(nanowire  +  H)  -  {F(nanowire)  +  £(H)}.  (4) 


The  results  of  A E  are  shown  in  Table  1.  In  these  results,  the  stabi¬ 
lization  energy  is  the  lowest  on  site  1  for  all  models.  Therefore,  a  H 
atom  is  adsorbed  on  the  top  of  B  (Al)  atoms.  In  the  following,  we 
only  consider  the  hydrogen  adsorption  on  the  top  of  Al  or  B  of  the 
ring  layer. 

Next  we  consider  the  adsorption  of  H2  molecules.  The  number 
of  H  atoms  is  even,  and  therefore  the  following  definition  of  the 
adsorption  energy  is  used: 


.  E(nanowire  +  N  x  H2)  -  {F(nanowire)  +  N  x  £(H2)} 

AE=  - iu - >  (5) 


N 
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Fig.  16.  The  largest  eigenvalues  of  stress  tensor  density  of  nanowire  models  after  hydrogen  adsorption.  Red  (blue)  region  means  positive  (negative)  stress  density.  The 
eigenvectors  are  normalized.  The  largest  eigenvaluesand  eigenvectors  of  stress  tensor  density  in  (a)  and  (b)  A1B  pentagonal,  (c)  and  (d)  A1B  hexagonal,  and  (e)  and  (f)  Al 
pentagonal  nanowire  models.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Table  2 

Geometry  parameters  and  the  stabilization  energy  of  H2  molecule  adsorption 


where  N  is  the  number  of  H2  in  the  system.  The  results  are  shown  in 
Table  2.  A E  is  calculated  on  the  condition  that  hydrogen  is  adsorbed 
on  all  adsorption  sites  as  shown  in  Fig.  11.  The  radius  of  B  ring,  rlt 
and  the  distance  between  B  and  H,  r2,  are  investigated  as  geometric 
parameters  with  keeping  their  symmetry  of  nanowires  (Fig.  11 ).  In 
the  case  of  the  pentagonal  A1B  nanowire  and  H  atoms,  n  and  r2 
are  determined  at  1.47  A  and  1.20  A,  respectively.  In  this  structure, 
the  bond  length  between  B  atoms  on  the  same  ring  is  1.73  A  and 
that  between  B  and  Al  atoms  is  2.50  A.  The  stabilization  energy 
A E  is  -2.14 eV.  In  the  hexagonal  A1B  nanowire,  r\  (r2)  are  1.68 
(1.20)  A,  the  length  between  B  (B  and  Al)  atoms  is  1.68  (2.34)  A, 
and  the  stabilization  energy  A E  is  -2.10eV.  In  both  models,  the 
bond  lengths  between  B  atoms  on  the  same  ring  are  longer  than 
those  of  A1B  nanowire  without  H  atoms.  The  stabilization  energy  is 
below  -2.0  eV.  On  the  other  hand,  in  the  case  of  Al  nanowire,  the 
bond  length  between  Al  atoms  on  the  same  ring  is  2.76  A,  which  is 
shorter  than  Al  nanowire  without  H  atoms.  The  stabilization  energy 
is  -0.5  eV.  Therefore,  A1B  nanowires  are  also  more  stable  than  Al 
nanowire  for  hydrogen  adsorption.  This  difference  of  the  hydrogen 
adsorption  energy  comes  from  the  effect  of  the  stabilization  of  the 
localized  electron  in  FIOMO.  Since  the  adsorption  energy  of  a  FI2 
molecule  adsorption  is  -1.83  eV  in  the  A1B2  bulk  surface  model  1 
and  -1 .12  eV  in  the  surface  model  2  in  the  case  of  A1B2,  these  results 
mean  that  A1B  nanowires  are  more  stable  than  A1B2  bulk  surface  for 
hydrogen  adsorption. 

3.2.2.  Electronic  structure 

The  DOS  and  band  structure  of  A1B  and  Al  simple  nanowires  after 
H  adsorption  are  shown  in  Fig.  12.  The  DOS  of  A1B2  bulk  surface 
models  after  hydrogen  adsorption  are  shown  in  Fig.  13.  The  DOS 
and  band  structure  of  all  models  are  metallic.  In  both  A1B  simple 
nanowires,  the  state  on  the  Fermi  level  consists  of  narrow  band. 
Thus  the  conductivity  of  A1B  simple  nanowire  is  different  from  the 
A1B2  bulk  and  surface. 

The  electron  density  on  the  xz-plane  and  xy-plane  are  shown 
in  Fig.  14.  Fig.  14(a),  (b),  (d),  (e),  (g)  and  (h)  are  the  sum  of  the 
occupied  states  and  Fig.  14(c),  (f)  and  (i)  are  the  highest  occupied 
state.  Fig.  14(a),  (d)  and  (g)  show  electron  density  on  xy-plane  and 
the  other  figures  show  that  of  xz-plane.  In  Fig.  15,  the  kinetic  energy 
density  and  the  tension  density  is  shown.  Fig.  15(a),  (b),  (d),  (e),  (g) 
and  (h)  are  the  sum  of  the  occupied  states  and  Fig.  15(c),  (f)  and  (i) 
are  the  highest  occupied  state.  The  tension  density  is  normalized. 
Fig.  15(a),  (d)  and  (g)  show  the  kinetic  energy  density  and  tension 
density  on  xy-plane  and  the  other  figures  show  those  of  xz-plane. 
Fig.  16  shows  the  largest  eigenvalue  and  eigenvector  of  stress  tensor 
density.  The  eigenvector  is  normalized.  Fig.  16(a),  (c)  and  (e)  show 
stress  tensor  density  on  xy-plane  and  the  other  figures  show  those 
of  xz-plane. 

First,  we  discuss  the  pentagonal  A1B  nanowire.  After  the  hydro¬ 
gen  adsorption  on  the  pentagonal  A1B  nanowire,  the  sum  of  the 
electron  density  in  Fig.  14(a)  and  (b)  appear  between  B  and  H.  In  the 
highest  occupied  state  (Fig.  14(c)),  the  electron  density  also  appears 
around  H  atoms.  In  the  kinetic  energy  density  (Fig.  15(a)-(c)),  the 
sum  of  the  all  occupied  states  shows  the  same  tendency  in  the 
electron  density.  The  kinetic  energy  density  regions  with  nega¬ 


tive  kinetic  energy  density  for  the  highest  occupied  state  appear  at 
the  center  of  the  five-membered  ring.  This  means  that  the  highest 
occupied  state  consists  of  FI  states.  Thus  the  peak  around  the  Fermi 
level  in  the  DOS  figure  is  the  electron  related  with  B-H  bonds.  In 
Fig.  16(a)  and  (b),  the  positive  region  and  the  spindle  structure  are 
seen  between  B  and  FI  and  between  B  atoms  on  the  same  ring.  This 
means  that  the  both  bonds  are  the  covalent  bond.  The  A1B  hexag¬ 
onal  nanowire  similar  properties  to  the  pentagonal  model  as  seen 
in  Figs.  12  and  14-16. 

In  Al  nanowire,  it  is  seen  in  Fig.  12  the  DOS  of  Al  nanowire  on 
the  Fermi  level  is  changed  by  the  hydrogen  adsorption.  The  elec¬ 
tron  density  of  the  sum  of  the  all  occupied  state  is  high  around 
FI  atoms.  In  the  highest  occupied  state,  the  electron  density  and 
the  kinetic  energy  density  are  localized  between  Al  and  H  and 
between  atoms  on  the  axis.  This  is  different  from  the  electron 
density  before  the  hydrogen  adsorption.  In  Fig.  16(e)  and  (f),  the 
positive  region  and  the  spindle  structure  are  seen  between  Al 
and  FI  and  the  negative  region  is  seen  between  Al  atoms.  This 
means  that  the  Al-H  bond  is  the  covalent  bond  and  Al-Al  bond  is 
metallic. 

4.  Conclusion 

We  performed  first-principles  calculations  for  A1B  nanowires 
to  study  the  transitions  of  their  electronic  structures  caused  by  FI 
atoms  adsorption.  The  difference  between  A1B  nanowires  and  Al 
nanowire  is  clearly  seen  in  the  electron  density  and  energy  density. 
The  electron  of  the  highest  occupied  state  appears  on  the  direction 
perpendicular  to  the  wire  axis  like  dangling  bonds.  Therefore,  it 
can  be  concluded  that  the  A1B  nanowires  do  not  have  the  conduc¬ 
tivity,  while  A1B2  bulk  has  it.  It  can  be  also  seen  that  A1B  nanowires 
are  more  stable  in  the  hydrogen  adsorption.  After  the  hydrogen 
adsorption  occurred,  H  atoms  make  the  covalent  bond  with  B  and 
Al.  There  is  no  significant  difference  of  the  electron  density  between 
the  hexagonal  A1B  nanowire  and  the  pentagonal  within  this  work. 
The  DOS  of  the  Al  nanowire  on  the  Fermi  level  is  decreased  by  the  H 
adsorption.  We  have  confirmed  that  conductivity  of  A1B  nanowires 
are  unchanged  with  hydrogen  adsorption  as  seen  from  the  electron 
density  and  the  energy  density.  The  hydrogen  adsorption  energy  of 
A1B  nanowire  models  is  larger  than  that  of  Al  nanowire.  In  this  point, 
A1B  nanowire  is  more  stable  in  the  hydrogen  adsorption. 
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